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Abstract 

The  assumption  of  local  thermal  equilibrium  (LTE)  between  the  gas  and  solid  matrix  within  porous  electrodes  of  solid  oxide  fuel  cells  (SOFCs) 
is  commonly  used  in  predictive  state-of-the-art  thermal- fluidic  models,  yet  its  validity  has  never  been  properly  established.  In  this  paper,  a  two- 
equation,  thermal  non-equilibrium  model  of  heat  transfer  between  the  gas  and  solid  phase  and  its  simple  scaling  analysis  are  used  to  estimate 
the  magnitude  of  the  temperature  difference  that  would  be  expected  between  the  gas  and  solid  phase,  and  thus  the  significance  of  local  thermal 
non-equilibrium  (LTNE)  effects.  A  simple  criterion  is  developed  for  determining  under  what  conditions  LTNE  effects  may  safely  be  neglected, 
and  when  they  are  likely  to  become  significant. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Thermal  modeling  is  a  vital  tool  in  the  design  and  develop¬ 
ment  of  SOFCs  due  to  the  close  coupling  between  temperature, 
reaction  kinetics,  ionic  resistance  of  the  electrolyte,  current  den¬ 
sity,  operating  voltage,  fuel  utilization,  and  flow  fields.  Accurate 
prediction  of  the  temperature  fields  within  SOFCs  is  not  only 
essential  to  predicting  and  optimizing  overall  cell  performance, 
but  is  also  a  prerequisite  to  mitigating  thermo-mechanical  degra¬ 
dation  and  failure  of  the  delicate  anode-electrolyte-cathode 
structure.  These  concerns  have  resulted  in  the  development  of 
highly  detailed  CFD/FEM-based  models  whose  ability  to  pre¬ 
dict  global  quantities  such  as  average  current  density,  voltage, 
and  temperature  has  been  well  established  through  experimen¬ 
tal  validation.  However,  lack  of  detailed  experimental  data  such 
as  local  temperature  fields  or  current  density  distribution  has 
prohibited  a  more  thorough  validation  of  the  numerical  models. 
Thus,  simplifying  assumptions  of  the  models  can  and  should  be 
challenged  from  a  physical  standpoint  in  order  to  maintain  rigor 
and  to  provide  the  required  careful  justification  (and  limits  of 
validity)  of  those  assumptions. 
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One  such  assumption  is  that  the  temperatures  of  the  gas  phase 
and  solid  matrix  within  porous  electrodes  are  locally  the  same, 

1. e.  the  gas  and  solid  are  in  local  thermal  equilibrium  (LTE).  A 
survey  of  the  relevant  literature  indicates  that  most  (possibly  all) 
of  the  current  models  make  this  assumption  with  little  or  no  jus¬ 
tification  given.  (Recently,  however,  Hwang  [1]  has  addressed 
this  issue  in  modeling  low  temperature  proton  exchange  mem¬ 
brane  fuel  cells.)  Three  typical  conditions  found  in  the  porous 
electrodes  of  SOFCs  bring  the  assumption  of  local  thermal 
equilibrium  into  question:  (1)  the  presence  of  volumetric  heat 
generation  in  the  medium,  (2)  very  low  Reynolds  number  flow, 
and  (3)  large  difference  in  thermal  conductivities  of  the  solid 
and  gas  phases.  Under  these  conditions,  the  gas  and  solid  tem¬ 
peratures  can  differ  greatly  and  local  thermal  non-equilibrium 
(LTNE)  modeling  is  often  required  [2].  To  carefully  address 
these  issues,  LTNE  effects  are  investigated  here  through  an  order 
of  magnitude  analysis,  which  results  in  the  development  of  a  sin¬ 
gle  LTNE  parameter  that  can  be  used  to  assess  the  validity  of 
the  local  thermal  equilibrium  assumption  in  SOFCs. 

2.  Model  formulation 

During  steady  state  operation,  all  of  the  heat  generated  within 
the  anode-electrolyte-cathode  structure  of  a  well-insulated  cell 
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Fig.  1.  Schematic  representation  of  the  TPB  region  of  the  porous  anode  and  solid  electrolyte  of  a  SOFC.  Oxygen  ions  migrate  from  the  cathode  side,  through  the 
electrolyte,  and  react  with  hydrogen  ions  adsorbed  on  the  surfaces  of  the  anode  to  form  water  vapor.  Heat  generated  due  to  electrochemical  reaction  and  ohmic  losses 
is  dissipated  by  three  mechanisms:  convection  to  the  gas  phase,  conduction  through  the  solid,  and  intra-  and  inter-pore  radiative  exchange. 


is  transferred  to  the  fuel  and  oxidizer  flow  streams,  which  carry 
it  out  of  the  cell.  In  the  porous  electrodes,  most  of  the  heat  is 
generated  near  the  electrode-electrolyte  interface  (see  Fig.  1) 
and  is  dissipated  by  (i)  conduction  through  the  solid  matrix,  (ii) 
heat  transfer  from  the  solid  to  the  gas  phase  by  convection  within 
the  porous  electrodes,  and  (iii)  advection  of  the  gas  through  the 
micro-pores  of  the  electrode  into  the  flow  channel.  These  pro¬ 
cesses  are  accounted  for  in  the  following  two-equation,  thermal 
non-equilibrium  model,  derived  from  conservation  of  energy  in 
the  gas  and  solid  phases,  and  applied  to  a  two-phase  representa¬ 
tive  elementary  volume  of  the  porous  electrode  [2] : 

0  —  V  •  (&s,eff  VTS)  +  hsfas(Tf  Ts ) 

+  (2 gen,  (solid  phase), 

V  •  ( pVcpTf )  =  V  •  (fef.eff  V7f)  -  hs{as(T{  -  Ts)  (gas  phase) 

(1) 

where  hs f  is  the  solid  to  gas  phase  convective  heat  transfer  coef¬ 
ficient,  as  the  specific  (per  unit  volume)  surface  area  within  a 
representative  volume  of  the  media,  p  the  density,  cp  the  specific 
heat,  and  V  the  velocity  of  the  gas  phase,  and  2  gen;  represents 
sources  of  volumetric  heat  generation  within  the  solid  phase.  The 
simplest  models  for  isotropic  effective  thermal  conductivity  of 
the  solid  and  gas  phases  are 

^s,eff  —  (1  o)ks,  kfeff  =  £kf  (2) 

respectively,  where  e  is  porosity.  This  space-homogenized, 
volume- averaging  model  is  most  commonly  used  in  modeling 
fuel  cell  porous  electrodes,  in  contrast  to  the  discrete-phase, 
microstructure-based  models,  in  which  two  phases  (gas  and 
solid)  are  spatially  separated  and  usually  analyzed  using  com¬ 
putationally  intensive  Lattice-Boltzmann  formalism. 

The  magnitude  of  the  difference  in  local  temperatures 
between  the  gas  and  solid  phases,  AT=  Ts  —  7f,  is  an  indicator 
of  the  validity  of  the  assumption  of  local  thermal  equilibrium 
(LTE).  Since  heat  generation  is  the  driving  force  behind  any 


thermal  non-equilibrium  in  SOFCs,  an  order  of  magnitude  anal¬ 
ysis  of  the  energy  conservation  equation  for  the  solid  phase  in 
Eqs.  (1)  yields  a  proper  estimate  of  the  magnitude  of  AT.  As 
a  first  approximation,  temperature  gradients  in  the  solid  (first 
term  on  the  right-hand  side)  can  be  assumed  small  compared  to 
the  other  two  terms  in  the  equation  and  can  be  neglected.  This 
requires  that  the  remaining  two  terms  balance  each  other,  and 
the  expected  temperature  difference  is  then  scaled  as 

O'" 

A  T-SJEL  (3) 

hsf(2s 

This  temperature  difference  is  the  greatest  possible  (upper 
bound)  that  can  be  expected  for  the  problem  at  hand,  because  it 
assumes  that  all  the  heat  generated  in  the  solid  is  dissipated  by 
convection  to  the  gas  phase.  For  verification,  examine  the  effects 
of  the  assumption  that  the  heat  conduction  term  is  negligible. 
As  solid  conduction  effects  become  significant,  this  provides 
an  additional  route  for  heat  transfer,  leading  to  a  reduction  in 
convective  heat  transfer,  and  thus,  reducing  the  temperature  dif¬ 
ference  between  solid  and  gas.  The  radiative  heat  transfer  is 
also  neglected  in  this  analysis.  In  similar  manner  to  heat  con¬ 
duction,  this  additional  heat  transfer  mechanism,  if  significant, 
would  provide  an  additional  route  for  heat  removal  via  intra-pore 
radiative  heat  exchange  within  the  solid  matrix  of  the  porous 
electrode,  thus  only  diminishing  the  expected  temperature  dif¬ 
ference  between  the  gas  and  solid  phases.  Therefore,  Eq.  (3)  is  a 
conservative  indicator  of  whether  or  not  LTNE  effects  should  be 
neglected.  The  remaining  task  is  to  estimate  the  magnitude  of  the 
heat  generation  term  and  the  interstitial  heat  transfer  coefficient, 
hv  =  hsfas ,  that  would  be  expected  in  the  SOFC. 

3.  Interstitial  heat  transfer  coefficient 

The  solid  to  gas  phase  heat  transfer  coefficient,  hs f,  depends 
on  the  thermal  conductivity  of  the  fluid,  the  Reynolds  number  of 
the  flow,  and  the  geometric  characteristics  of  the  medium.  The 
specific  surface  area,  as ,  also  depends  on  the  structure  of  the 
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solid  matrix.  As  a  first  approximation  for  this  order  of  magni¬ 
tude  analysis,  the  structure  of  the  porous  electrode  material  can 
be  assumed  isotropic,  consisting  of  roughly  spherical  particles 
of  diameter,  dv.  This  is  a  conservative  approximation — more  so 
than  considering  a  continuous  solid  phase,  which  would  enhance 
conductive  heat  transfer,  thereby  reducing  LTNE  effects.  The 
SOFC  electrodes  have  pores  or  voids  of  diameter,  dv,  which  has 
been  experimentally  measured  as  0. 3-0.4  p,m  for  anode  materi¬ 
als  [3] .  (For  simplicity,  we  assume  that  pore  size  in  the  cathode  is 
of  the  same  magnitude.)  The  relationship  between  specific  sur¬ 
face  area  and  void  diameter  given  by  Vafai  [4]  is  as  =  As/dw.  This 
yields  a  specific  surface  area  of  the  order  of  4  x  106m  ^Parti¬ 
cle  diameter  can  be  estimated  from  the  equation  for  the  specific 
surface  area  for  either  spherical  particles  [5]: 


as 


6(1  -  e) 

dp 


or  for  sintered  metal  beads  [6] : 


as 


20.346(1  -  s)s 2 

dp 


yielding  0.9  and  1.1  |jim,  respectively.  Although  the  exact  value 
of  these  parameters  is  uncertain  for  the  available  porous  elec¬ 
trode  materials,  it  is  clear  that  as  is  proportional  to  1  /dp,  and 
particle  and  void  diameters  are  the  same  order  of  magnitude. 

Another  area  of  uncertainty  in  determining  interstitial  heat 
transfer  coefficient  is  the  Reynolds  number,  Red  =  pfudvliif,  of 
the  gas  flow  through  the  micro-pores  of  the  electrodes.  In  gen¬ 
eral,  it  is  reasonable  to  expect  that  the  transport  through  the 
electrodes,  especially  near  the  electrode-electrolyte  interface, 
is  dominated  by  diffusion  [3,7].  In  this  region,  known  as  the 
triple  phase  boundary  (TPB)  region,  the  electrochemical  reac¬ 
tion  takes  place  at  the  junction  of  catalyst,  electrolyte  and  gas 
phase.  As  a  first  estimate,  the  order  of  magnitude  of  the  aver¬ 
age  velocity  here  can  be  found  from  conservation  of  mass  and 
Faraday’s  law,  which  relates  moles  of  reactants  or  products  to 
electrical  current.  Specifically,  for  a  given  current  density,  the 
required  mole  flux  of  the  relevant  species  (e.g.  H2)  to  the  TPB 
is  given  by,  h  =  i/2F,  where  i  is  the  current  density  of  the  cell 
(A  m-2),  F  is  Faraday’s  constant  (96,500  C  mol-1  e-1),  and  the 
2  appears  because  2  electrons  are  produced  by  the  balanced 
electrochemical  reaction.  For  a  representative  current  density  of 
104  A  m-2  the  corresponding  mass  flux  of  H2  to  the  TPB  region 
is  of  the  order  0.1  gm-2  s_1.  Dividing  by  density  and  porosity 
yields  a  gas  phase  velocity  on  the  order  of  0.01  ms-1.  Given 
such  small  velocity  and  using  the  previously  estimated  particle 
diameter  and  appropriate  fluid  properties,  it  becomes  apparent 
that  Red  1  in  the  electrochemically  active  region  of  the  porous 
electrode. 

A  very  common  correlation  for  hs f ,  developed  by  Wakao  and 
Kaguei  [8]  from  an  extensive  collection  of  experimental  data 
(for  spherical  particles)  is 


2+1.1  Re°d6Prf33 


where  kf  is  the  thermal  conductivity  of  the  fluid  and  Pr f  is 
the  Prandtl  number  of  the  fluid.  This  correlation  is  valid  over 
a  large  range  of  Re d,  although  little  data  exists  for  very  low 
Reynolds  numbers.  However,  the  correlation  does  asymptoti¬ 
cally  approach  the  correct  theoretical  value  of  h$ f  in  the  limit  of 
Re d  — >  0  (pure  conduction  through  the  gas  phase).  (Other  fre¬ 
quently  used  correlations,  such  as  that  of  Kar  and  Dybbs  [9], 
do  not  appear  to  be  valid  at  this  limit  [2].)  Using  the  fact  that 
Re d  1  as  previously  discussed,  Eq.  (6)  is  reduced  to 

2kf 

ftsf  -  ~r  (?) 

dv 

This  is  the  most  conservative  estimate  of  hs f  since  it  considers 
only  heat  transfer  by  conduction  in  the  gas.  Any  increase  in 
Re d  would  only  increase  the  heat  transfer  coefficient  from  this 
baseline  value. 

Combining  Eq.  (7)  with  Eq.  (4)  or  (5)  yields  a  lower  bound 
estimate  of  the  interstitial  heat  transfer  coefficient,  hY.  A  typical 
value  of  this  parameter  for  air,  assuming  electrodes  with  1  jam 
particles  and  porosity  of  30%  is  on  the  order  of  101 1  W  m~3  K_1 . 
This  very  high  volumetric  heat  transfer  coefficient  is  due  to  effi¬ 
cient  heat  transfer  within  the  microscale  porous  structure  of  the 
electrodes  because  of  the  very  large  specific  surface  area.  Since 
the  calculations  are  sensitive  to  particle  diameter  and  porosity, 
results  will  be  presented  for  a  range  of  these  parameters. 

4.  Estimate  of  volumetric  heat  generation 

An  upper  limit  for  the  magnitude  of  the  volumetric  heat  gen¬ 
eration  term,  Q'"q n  in  Eq.  (3)  is  found  from  an  examination 
of  the  thermodynamic  effects  associated  with  electrochemical 
reactions  in  the  fuel  cell  and  the  ohmic  heating  effects  in  the 
ionic  conducting  electrolyte  material.  Although  the  processes 
by  which  heat  is  generated  are  complex  and  coupled,  the  overall 
energy  balance  of  the  cell  effectively  provides  an  estimate  of 
the  upper  limit  of  heat  generated  without  requiring  an  in-depth 
analysis  of  these  processes.  The  total  heat  generation  for  a  sin¬ 
gle  cell  can  be  calculated  based  on  a  given  current  density,  and 
that  entire  quantity  of  heat  is  then  applied  to  the  TPB  region 
(assumed  to  be  a  10  |xm  thick  layer  next  to  the  electrolyte  [3]) 
on  only  one  side  of  the  electrolyte  in  order  to  estimate  the  upper 
bound  of  volumetric  heat  generation. 

The  reversible  electrochemical  reaction  for  the  hydrogen- 
fueled  cell  [10]  is 

H2  +  \02  +  2e-  =  H20  +  2e-  (8) 

The  limit  of  the  energy  available  for  electrical  work  is  the  Gibbs 
free  energy,  given  at  standard  temperature  (298  K)  and  pressure 
(1  atm)  as 

AG°  =  A H°  -  T°AS°  =  —nFE°  (9) 

where  n  is  the  number  of  electrons  (2)  per  mole  of  reactant,  F 
is  Faraday’s  constant  and  the  °  indicates  standard  values.  For 
operating  conditions  other  than  standard,  the  Gibbs  free  energy 
is  given  by 


(6) 


AG  =  AH  —  TAS  ~  —nFEidezi 


(10) 
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where  £ideai  is  the  ideal  (reversible)  operating  voltage  of  the  cell 
at  operating  temperature  and  pressure.  Adding  Eqs.  (9)  and  (10) 
and  solving  for  TAS  yields  the  total  reversible  heat  generated 
(kJ  mol-1)  by  the  reaction: 

Q  rev  =  -TAS 

=  -T°AS°  +  nF(E°  -  £ideai)  +  (A H°  -  AH)  (11) 


Faraday’s  law  gives  the  rate  of  heat  generation,  and  this  is  written 
on  a  volumetric  basis  as 


Q 


rev 


L(1  -  e) 

T°AS°  +  (AH°  -  AH) 


x 


nF 


T  (E  F'ideal) 


(12) 


where  i  is  the  current  density,  L  the  length  scale  (thickness)  of 
the  heat  generating  region,  and  s  is  the  porosity. 

Because  the  cell  does  not  operate  ideally  (reversibly),  the 
operating  voltage  is  lower  than  the  ideal  voltage.  This  voltage 
drop  is  primarily  due  to  the  ionic  resistance  of  the  electrolyte, 
activation  overpotential,  and  concentration  overpotential,  result¬ 
ing  in  additional  (irreversible)  heat  generation  which  can  be 
calculated  as 


O'"  = 

V^irrev 


L(l-e) 


(-^ideal  ^oper) 


(13) 


Thus,  by  adding  Eqs.  (12)  and  (13),  the  total  volumetric  heat 
generation  is  expressed  as 


Q  total 


L(l-e) 

T°AS°  +  (AH°  -  AH) 


x 


nF 


+  (E  E0  per) 


(14) 


For  example,  consider  a  cell  operating  at  E0per  =  0.5  V  and  a  cur¬ 
rent  density  of  i  =  104  Am-2.  For  the  reaction  given  in  Eq.  (8), 
the  standard  potential,  E° ,  is  1.22  V,  T°AS°  is  — 13.2  kJ  mol-1 
(for  water  vapor  as  a  product),  A H°  is  —  241.8kJmol-1,  and 
AH  (evaluated  at  1000  K)  is  — 247.9  kJ  mol-1  [11].  If  all  the 
generated  heat  is  confined  to  the  10  p,m  thick  TPB  region  and 
assuming  a  porosity  of  0.3,  the  volumetric  heat  source  term  is  of 
the  order  109  W  m-3 .  Combining  this  result  with  the  estimate  for 
interstitial  heat  transfer  coefficient  found  in  Section  3,  Eq.  (3), 
gives  an  estimate  of  the  local  temperature  difference  between 
gas  and  solid  phases  to  be  of  the  order  10-2  K,  i.e.,  negligible. 


5.  Estimate  of  LTNE  effects 


Having  developed  expressions  for  heat  generation  and  inter¬ 
stitial  heat  transfer  coefficient,  the  LTNE  criteria  given  in  Eq. 
(3)  can  now  be  used  to  estimate  the  maximum  temperature  non¬ 
uniformity  between  the  solid  and  gas  phases.  However,  to  use 
Eq.  (14)  to  calculate  heat  generation,  the  operating  temperature, 
voltage  and  current  must  be  known.  If  these  are  not  known,  an 
approximate  upper  bound  is  found  by  setting  the  operating  volt¬ 
age  in  Eq.  (14)  to  zero  (i.e.,  cell  is  short-circuited  and  drawing 
maximum  current),  and  neglecting  the  first  term  in  the  brackets. 


Fig.  2.  The  maximum  expected  temperature  difference  between  the  solid  and  gas 
phases  of  the  porous  electrode  as  a  function  of  current  density  for  various  porosity 
and  particle  diameter  of  the  porous  medium  (computed  using  representative 
L  =  10  |xm  and  thermal  conductivity  of  air,  kf  =  0.5  W  m- 1  K- 1 ). 


This,  along  with  Eqs.  (4)  and  (7)  are  substituted  into  Eq.  (3), 
yielding  the  upper  bound  of  gas-to-solid  temperature  difference 
that  can  be  expected  in  SOFC  electrodes: 

id2 

ATmax-O-lx  —f  -  (15) 

kfL(l  -  sy 

This  parameter  allows  the  user  to  determine  if  LTNE  effects 
can  safely  be  ignored,  or  if  detailed  non-equilibrium  thermal 
modeling  is  required.  Fig.  2  shows  a  plot  of  ATmax  versus  elec¬ 
trode  particle  diameter  for  various  values  of  current  density  and 
porosity.  Although  our  assumptions  have  erred  on  the  side  of 
maximizing  the  expected  temperature  difference,  it  should  be 


0  2  4  6  8  10 

TPB  thickness,  L  [pm] 


Fig.  3.  The  maximum  expected  temperature  difference  between  the  solid  and 
gas  phases  of  the  porous  electrode  as  a  function  of  triple  phase  boundary  (TPB) 
thickness.  The  computations  are  performed  for  representative  parameters  of 
i=  1000  mA  cm-2,  c/p  =  1  pan,  £  =  0.3,  and  kf  =  0.5  W  m_l  K-1. 
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noted  that  local  current  density  could  be  significantly  higher 
than  the  average  value  used  in  Eq.  (15)  leading  to  local  “hot 
spots”  with  increased  LTNE  effects.  In  addition,  the  thickness 
of  the  TPB  region,  L,  may  be  much  smaller  than  10  [xm  used  in 
computing  results  presented  in  Fig.  2,  resulting  in  an  increased 
ATmax.  To  illustrate  the  magnitude  of  effect,  Fig.  3  is  a  plot 
of  ATmax  versus  L  for  typical  conditions  of  i=  1000  mA  cm-2, 
dv  =  l  pan  and  £  =  0.3.  Although  ATmax  does  indeed  increase 
with  a  decrease  in  L,  the  temperature  difference  is  still  insignifi¬ 
cant  (less  than  one  degree).  This  result  should  be  used  with  care 
and  only  as  a  qualitative  trend  indicator,  since  the  basic  validity 
of  volume-averaging  model  formulation  becomes  questionable 
as  the  TPB  region  size  L  becomes  comparable  or  smaller  than  the 
elementary  representative  volume  used  for  deriving  Eq.  (15).  It  is 
left  to  the  discretion  of  the  user  to  properly  interpret  the  results 
of  Eq.  (15)  and  use  relevant  local  values  of  /,  kf,  L,  and  £  in 
assessing  the  validity  of  the  LTE  assumption  for  specific  cases. 

6.  Concluding  remarks 

An  order-of-magnitude  analysis  of  the  energy  conservation 
equation  for  the  solid  phase  of  the  SOFC  porous  electrodes  is 
used  to  estimate  the  expected  temperature  difference  between 
the  gas  and  solid  phases.  This  requires  an  estimate  of  the  vol¬ 
umetric  heat  transfer  coefficient  between  the  two  phases,  and 
an  estimate  of  the  volumetric  heat  generation  within  the  solid 
phase.  Based  on  this  analysis,  a  simple  relationship  is  developed 
for  estimating  the  maximum  possible  temperature  difference 


that  exists  between  gas  and  solid  in  SOFC  porous  electrodes.  If 
the  expected  temperature  difference  (or  error)  is  small  enough 
to  be  tolerated,  then  the  assumption  of  LTE  is  justified.  How¬ 
ever,  caution  should  be  exercised,  as  this  analysis  is  specific 
to  consideration  of  the  idealized,  isotropic,  porous  electrodes 
of  hydrogen-fueled  SOFCs  and  numerous  simplifying  assump¬ 
tions  have  been  made  in  estimating  local  cell  parameters  and 
structural  characteristics. 
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